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We study two non-equilibrium work fluctuation theorems, the Crooks' theorem and the Jarzynski 
equality, for a test system coupled to a spatially extended heat reservoir whose degrees of freedom 
are explicitly modeled. The sufficient conditions for the validity of the theorems are discussed in 
detail and compared to the case of classical Hamiltonian dynamics. When the conditions are met 
the fluctuation theorems are shown to hold despite the fact that the immediate vicinity of the test 
system goes out of equilibrium during an irreversible process. We also study the effect of the coupling 
to the heat reservoir on the convergence of (exp(— to its theoretical mean value, where W is 
the work done on the test system and f} is the inverse temperature. It is shown that the larger the 
local heating, the slower the convergence. 
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I. INTRODUCTION 

The recently discovered non-equilibrium fluctuation 
theorems deal with large deviations and symmetries 
of quantities [iI_3L such as entropy production d, Q], 
heat exchange [inSl J or work performed 0, [1| during 
an irreversible process. Two examples which have re- 
ceived much attention are the Jarzynski equality and the 
Crooks' theorem. In the Jarzynski equality one considers 
a test system (TS) which is initially in equilibrium with 
a heat reservoir at temperature T. The test system is 
then driven to another state by varying an external pa- 
rameter A over a finite time interval < t < t using a 
fixed protocol X{t). In this process, a history dependent 
work, W, is performed. The Jarzynski equality [31 states 
that the average of e"*^^ over histories satisfies 
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where (3=1 /T is the temperature of the system when it 
is in equilibrium and AFts ~ Fts{B) — Fts{A) is the 
free energy difference between the equilibrium states, A 
and i3, of the test system at the initial and final values 
of A, respectively. 

The Crooks' theorem is a statement about the prob- 
ability distribution of work performed in an irreversible 
process Q . As for the Jarzynski equality, one considers a 
test system in contact with a heat reservoir at tempera- 
ture T. Initially in an equilibrium state A with A = A(0), 
the test system is driven, during the time interval (0,r), 
by the protocol \{t) (the 'forward process'). During this 
process, varying amounts of work may be performed, de- 
pending on the trajectory the system takes in state-space 
as A is varied; the likelihood that work W is performed is 
specified by a probability distribution Pf{W). One then 
considers the reverse process in which the system starts 
in the equilibrium state B with A = A(t), and is driven 
by the time-reversed protocol \{t) = A(t — t) in the time 
interval (0, r). In this process the probability of work W 
being done is specified by the distribution PuiW). The 



Crooks' theorem states that 
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Note that the Crooks' theorem is a stronger statement 
than the Jarzynski equality; indeed the latter can be ob- 
tained by integrating the former 

As has been discussed extensively in the past 0, H, 
[lo| , both relations rely on rare events in which the work 
performed deviates significantly from the average. Thus, 
the relevance of Eqs. ([I]) or as a tool for measuring 
free-energy differences is restricted to small systems, for 
which the probability of sampling atypical events is non- 
negligible even for experimentally reasonable realizations 
of the protocol [ll| . 

Since their discovery a decade ago, the fluctuation the- 
orems have been proven for both stochastic [1, H, 
and Hamiltonian dynamics [3, [ij, |T^. Although nu- 
merous proofs have been given, the precise assumptions 
neces sary for their validity has been the subject of de- 
bate [la [IE [III- In large measure, this lack of clarity 
revolves around the importance of local heating, to wit, 
that near the region of contact with the test system the 
heat reservoir is ziseZ/ driven out of equilibrium during the 
experiment. For example, proofs relying on Markovian 
dynamics [1| assume that the test system obeys detailed 
balance with an ideal heat reservoir, that is, the tran- 
sition rate, wjf from conflgurations i to j of the test 
system satisfies 
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where Ei and Ej are energies of the configurations i and 
j of the test system, respectively. Similarly, for Langcvin 
dynamics one typically assumes white noise satisfying the 
fluctuation-dissipation relation [I^ . Thus, these dynam- 
ics implicitly assume that the heat reservoir is always in 
thermal equilibrium [23|. While local heating may, in 
practice, be small, it is nonetheless a fundamental issue 
which goes to the heart of the validity of the fluctuation 
theorems 
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The aim of this paper is to elucidate the requirements 
for the validity of the Crooks' theorem and the Jarzynski 
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equality in the presence of local heating. While several 
work have considered a restricted class of non-Markovian 
sources of noise acting as the heat reservoir [l^l , we take 
a different approach. To this end, we consider the test 
system, on which the driving is performed, to be a por- 
tion of a larger isolated "combined system" consisting 
of the test system + heat reservoir, with no assumption 
made about the heat reservoir degrees of freedom. We 
first review results for classical Hamiltonian dynamics 
and their limitations, and then present a new proof for 
Markovian stochastic dynamics which requires less strin- 
gent restrictions than their classical Hamiltonian coun- 
terparts - specifically for classical Hamiltonian dynam- 
ics the motion in phase space is incompressible while for 
stochastic dynamics there is no such restriction. We show 
that in order for Crooks' theorem to hold, it is sufficient 
that the following three conditions are satisfied: 

1. The closed "combined system" (CS), consisting 
of the test system -I- heat reservoir, obeys time- 
reversal symmetric dynamics when the external pa- 
rameters (such as A) are constant. 

2. When in thermal equilibrium, extensive properties 
of the systems are additive, e.g., the free energy of 
the CS is the sum of the free energies of the test 
system and the heat reservoir: Fes = Fts + Fb- 

3. The control parameters (such as A) directly couple 
only to the test system (though their influence may 
surely be felt in the entire CS). 

We note that when condition (2) docs not hold, modified 
relations may still be obtained; we comment on this in 
the body of the paper. Since we make no assumptions 
about the state of the heat reservoir, this demonstrates 
that local heating invalidates neither the Jarzynski equal- 
ity nor the Crooks' theorem. This notwithstanding, local 
heating strongly impacts on the statistics needed to ac- 
curately estimate free energy differences, and we address 
this issue numerically. 



II. HAMILTONIAN DYNAMICS 

In the following we present a proof of the Crooks' theo- 
rem for an isolated system obeying classical Hamiltonian 
dynamics, along the lines of Reference This will al- 
low us to emphasize the differences between stochastic 
and classical Hamiltonian dynamics. We wish to stress 
that for closed systems (such as the CS above) obeying 
classical Hamiltonian dynamics, the Crooks' theorem re- 
lies on the incompressibility of trajectories in phase space 
(Liouville's theorem) and microscopic time reversibility. 

Consider the allowed phase-space manifolds T,e and 
^E+w of the CS with energies E and E + W respectively 
(see Figll]). As the system is driven using the forward 
protocol A(i), a subset cr^ of is mapped into a subset 
te+w of S E+w - for these trajectories an amount of work 
W is performed. Other subsets of arc mapped by the 
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FIG. 1: Hamiltonian dynamics: A schematic illustration of 
the mapping of sets of states, initially with energy E, un- 
der the protocol \{t). (a) In an irreversible process, different 
amounts of work {e.g., W, W', and W") may be performed, 
meaning that states from the manifold E_b may be mapped 
to manifolds of different energies. Note that due to the deter- 
ministic dynamics as and a'^ have no common elements and 
depend explicitly on the protocol X{t). (b) For a reversible 
process, the manifold E_b is mapped into the manifold T^e+w 
where volumes of both the manifolds are same. 



dynamics to other manifolds of different energies, for ex- 
ample 'Se+w'j as indicated in the figure. Denote the 
phase space volumes of aE and aE+w by uje and uje+w 
and those of and T^e+w and fls+Wi respec- 

tively. The incompressibility of phase space trajectories 
guaranties that uje = ^^e+w- Microscopic time reversal 
symmetry (Condition 1) implies that if the driving proto- 
col is reversed [l^, then cTe+w ^^^^ precisely mapped 
onto cr^, where (t'^e+w ^^'^ '^e ^^'^ time- reversed im- 
ages (i.e., all momenta reversed) of (Je+w and ct^;, re- 
spectively. These clearly all have the same volume in 
phase space. In the special case when the protocol is a 
reversible process, the entire manifold "Se will be mapped 
onto the entire manifold "Ee+w, where W is the energy 
added to the system in the reversible process, with no 
change in entropy. 

Since in equilibrium all states on a constant energy 
manifold are equally probable, the probability of per- 
forming work W in the forward protocol is given by 
Pf{W) = oje/^e- Similarly, for the reverse process 
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Pr{—W) = oje+w/^e+w- We expect that in the hmit 
of a large heat reservoir > PF, PpiW) and Pr(-PF) 
will not depend on E, which, together with loe = ^e+w 
gives 

Pf{W)/Pr{-W) = Qe+w/^e 

= gScs{E+W)-ScsiE) _ (-^-j 

Here we have identified il,E+w = e^'^-^'-^"'"^^ (and like- 
wise for ^e), where Scs{E + W) is the entropy of the 
CS with energy E + W. We may now couch this result in 
the standard form (Eq. ^ : First expand to first order in 
W and use AFcs = W- TAScs- Next, use Condition 
2 that Fes = Fts + Fb, and Condition 3 that Fb is 
independent of A (and so AFcs = AFts)- 

Note that for classical systems the condition Fes = 
Fts + Fb is rather strong. It implies that either a ther- 
modynamic limit of the test system is taken or that the 
interaction between the test system and the heat reser- 
voir is weak. The latter would imply negligible local heat- 
ing while the former would restrict the result to large test 
systems. If neither of these cases holds, we can not de- 
fine the test system as a distinct object. In this case, as 
evident from the above proof, the Crooks' theorem still 
holds with AFts replaced by AFcs implying that one 
has to consider the free-energy difference of the full com- 
bined system. Its value is expected to be dependent on 
the interaction energy between the test system and the 
heat reservoir. As we argue below, for stochastic systems 
these restrictions are less severe and one can have strong 
local heating even when Fes = Fts + Fb ■ 



III. STOCHASTIC DYNAMICS 
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FIG. 2: Stochastic dynamics: A schematic illustration of the 
mapping of sets of states, initially with energy E, under the 
protocol X{t). In this case, the dynamics do not map a given 
initial to a single final state. Time-reversibility only implies 
that a forward trajectory has a corresponding time-reversed 
trajectory of equal statistical weight, (a) In an irreversible 
process, different amounts of work {e.g., W, W' , and W") 
may be performed, meaning that states from the manifold E_b 
may be mapped to manifolds of different energies. In this case, 
(T_B and (T^ may have common elements, (b) For a reversible 
process, the manifold E_b is mapped into the manifold T.e+w. 
This notwithstanding, a given initial state may evolve along 
many different trajectories Y{t), as indicated. 



Stochastic dynamics arise either due to incomplete 
knowledge of some degrees of freedom of the system and 
the heat reservoir or due to quantum effects (or both). 
Although the choice of an effective stochastic dynamics 
may be justified from the underlying microscopic classical 
or quantum dynamics [U, [H, [2^ , it is not obvious that 
the fluctuation theorems are valid when only the effective 
stochastic dynamics is provided. For example, noise of 
quantum origin (or in a more general purely mathemat- 
ical settings) is not captured by Hamiltonian dynamics 
and thus invalidates the proof for the Hamiltonian case. 
The central difference between the stochastic case and 
that of classical Hamiltonian dynamics is that here the 
evolution of a volume in phase space is not incompress- 
ible. This means that a given initial state may evolve 
under the dynamics into many final states. Employing 
the same notation as in section II, this difference is illus- 
trated in Figure [21 which shows that the volume oje '^s 
not, in general, equal to oje+w- 



A. Crooks' Theorem 

We now turn to the proof of the Crooks' theorem for 
a stochastic time reversible dynamics. To understand 
irreversible and dissipative physics in a stochastic set- 
tings, there are various ways one can model a system 
[21| . We do in the following fashion. We consider an iso- 
lated CS comprised of a test system coupled to a reser- 
voir. Per Condition 1, we assume that the CS evolves 
through Markovian dynamics obeying microscopic time 
reversibility. Formally, this means that the transition 
rate Wj,i {wi_j) from configuration i to j (j to i) of the 
CS satisfies 

Wj^i = Wij ; (5) 

This is just the statement of detailed balance for an iso- 
lated system [^J] . We stress that this docs not imply Eq. 
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|3l that is, the effective heat reservoir is neither Marko- 
vian nor in quasi-equilibrium, and the above condition is 
a much weaker assumption than that made in References 

Let us now consider the 'forward' process in the time 
range {—oo < t < oo), where A(t) = Xa for (— oo,0), 
Xb for (r, cxd), and varies in a fixed fashion in (0, r). At 
t — ~oo the CS is in an cquiUbrium state A, and at i = oo 
the CS comes to an equihbrium state B. For a given pro- 
tocol X{t) we indicate the probabihty of a given trajectory 
Y{t) of the CS in the forward process by the functional 
VF[Y{t),X{t)]. Y{t) = {x{t),X{t)} includes aU the de- 
grees of freedom of the CS, that is those of the test system 
{x{t)) and of the heat reservoir {X{t)). The hmctional 
VR[Y{t), X{t)] denotes the probability of the correspond- 
ing time-reversed trajectory Y{t) = Y{—t) and time re- 
versed forcing protocol, X{t) = X{—t), in which the CS 
begins in state B at t = —oo and ends in ^ at i = cxd. 
Eq. [5] implies that 

VF[Y{t),X{t)]^VR[Y{t),X{t)] (6) 

for any trajectory. This simply says that the probabil- 
ity of any trajectory in the configuration space of the 
combined system is equal to the probability of its' cor- 
responding time reversed trajectory. Work performed is 
defined as the energy difference between the two ends 
of the trajectory, i.e., [Hcs{Yf, A/) - Hcs{Yt, AJ] where 
Hcs{Y{t), X{t)) is the total energy of the combined sys- 
tem at any time i, subscript i and / denote the initial and 
final points of the trajectory. We denote the functional 
WF[Y{t), X{t)] as the work performed along the for- 
ward trajectory, and WTi[Y{t), X{t)] is that for the corre- 
sponding time-reversed trajectory, with WF[Y{t), X{t)] — 
—yVji\Y{t),X{t)]. The probability that an amount of 
work W is performed in the forward and reverse pro- 
tocols, respectively, are given by 

PpiW) = J P^^Y_oo)Vf[Y, X]diWF - W)VY (7) 

and 

PRiW) = J P^''iYoo)VR[Y,X]SiWR - W)VY (8) 

where T>Y = Il^-oo "^^i denotes the integration over 
all trajectories in the configuration space, P'^{Y-ao) and 
P'^iYrxi) are the initial and final equilibrium distribu- 
tions of the CS, and we have omitted the explicit trajec- 
tory dependence of W_F[y, A] and >V/j[F, A]. 

The CS is in equilibrium for t = ±oo, and an amount 
of work W is performed only during the interval (0,r). 
For t = ±oo the CS is described by a microcanon- 
ical ensemble, and the probability of any microstate 
Y = {x, X} is the inverse of the total number of mi- 
crostates, e-'S'cs(*=±oo)^ expand Scs{E + W) to 

first order in W and use IS.Fcs = W — T^^Scs where the 
inverse temperature = dScs{E)/dE. According to 



Condition 3, A couples directly only to the test system, 
while Condition 2 implies that the free energy of the CS 
is the sum of the free energies of the reservoir and test 
system. Thus we have that 

PaO^-oo) _ P{Wf-AFts) (q\ 

Ps'iYoo) ' 

where Wj? is the work performed along the forward path 
and AFts = Fts{B) — Fts{A) is the free-energy differ- 
ence of the test system at < = oo and t = — oo, that is, 
for the values Xb and Xa [l^l ■ Putting Eq. [S] and Eq. [5] 
in Eq. [7l and using Wf = — W/j, wc get 

Pf{W) = j e-P^^''^^^^^^P^\Y^)VR[Y,X] 

x5{Wr + W)VY 
^^P{w-afts) j p^'1{y^)Vr[Y,X]5{Wr + W)VY 

= e'^'^^^^^-^'^PRi-W) (10) 

which is obtained after comparing with Eq. [S] This 
proves the Crooks' theorem for stochastic dynamics. 
Note that also here we assume that Pf{W) and Pr{W) 
are independent of E in the limit of large heat reser- 
voir. Finally, where Conditions 2 and 3 do not apply 
similar considerations to those discussed for Hamiltonian 
dynamics hold. 

Note that the conditions on the validity of Crooks' the- 
orem are weaker here than for Hamiltonian dynamics. 
For Hamiltonian dynamics, Liouville's theorem and Con- 
dition 1 imply time reversibility in the sense that each mi- 
crostate in the initial manifold maps under the action of 
A(t) onto a single final microstate, and under the reverse 
protocol the final microstate, with all momenta reversed, 
in turn maps back to the original microstate with all mo- 
menta reversed. In the stochastic case the constraint is 
weaker in the sense that an initial microstate may end up 
in any one of many final microstates; all that is required 
for the proof is that the weights of path connecting states 
in the forward and reverse direction are equal. 

In contrast to classical Hamiltonian dynamics, for 
stochastic dynamics there is no connection between lo- 
cal heating and the interactions between the test-system 
and the heat reservoir. Indeed, the energy function of 
the stochastic system may contain no interaction term 
between the test system and the heat reservoir. How- 
ever, this does not mean that there is no heat transfer, 
and consequently local heating can occur. This is illus- 
trated in the numerical example of Sec. HI. 

B. Jarzynski Equality 

Although the Jarzynski equality can be derived by di- 
rect integration of Crooks' Theorem, it is instructive to 
present a proof for stochastic dynamics in the spirit of 
the previous section. We discretize the process and con- 
sider a general protocol {Aq, Ai, . . . , A„^ }, where A^ is the 



value of A for the k update interval. The process con- 
sists of changing the value of A (while leaving Y = {x,X} 
unchanged), updating the CS, changing A again, etc. De- 
noting the values of the degrees of freedom of the CS at 
step k by Yk = {xk,Xk), this process is specified by the 
sequence of transitions {Yo,Xo) (yo,Ai) — > (Yi, Ai) 
(Yi, A2) . . . . The transition probability matrix between 
the states (Yfc_i, Afe) and (Yk,Xk), which has an implicit 
dependence on time through the changing parameter A, 
is denoted by Wk^k-i- Condition 1, with the Markovian 
dynamics for the entire CS implies that 

I Wk+i.kPcUy>^,\)dYk = P^'^siYk+uX) . (11) 

(Condition 1 guarantees that dynamics for a given value 
of A only connects microstatcs in the equilibrium ensem- 
ble.) 

Conditions 2 and 3 imply that the work at time step 
k is defined as AWk = H{xk, Xk+i) — H{xk,Xk) where 
H(x, A) is the energy of the test system in microstate x 
at parameter value A. The total work performed in the 
process is thus given hy W = ^^^fc- Condition 2 

implies that 

P^%{x, A) = cxp{~(3H{x, X))/Z , (12) 

where Z is the partition function of the test system. 
From Eqs. [TT] and fT2l we have 

j e-^^^^Wk+,.kP^UYk)dYk = ^P^%{Yk+,) , 

where Zk is defined in Eq. [T^with A = Afc. Using this 
repeatedly in the relation 

P^siYo,Xo) n e-P^'^^Wk+iM , 

fc=0 

we readily obtain Eq. [T]by identifying Fk = — In(Zfc), 
the free energy of the test system for A = A^. 

IV. NUMERICAL MODEL 

In the preceding sections we have shown that local 
heating does not invalidate either the Crooks' Theorem 
or the Jarzynski Equality. This notwithstanding, local 
heating has a crucial effect on the convergence of the re- 
sults. We illustrate these results using a simple stochastic 
model. We first show numerically the validity of both the 
Jarzynski equality and the Crooks' theorem for stochas- 
tic dynamics respecting Eq. O This allows us to study 
the effect of local heating on the typical number of ex- 
perimental runs needed to estimate the average in Eq. [1] 
It is shown that the larger the local heating, the greater 
the number of runs needed. 

Consider a one-dimensional [2^ lattice of L -|- 1 sites, 
with sites labeled i = 0,...,L; this is our CS. The test 
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FIG. 3: Schematic diagram of the CS. The test system occu- 
pies site i = 0, and the sites of the heat reservoir are labeled 
hy i = 1,2,3...L where L is very large. Interactions are 
between neighboring sites. 



system is placed at site i ~ 0, and the heat reservoir 
occupies sites 1 < i < L, as indicated in Fig. [31 At each 
site we define a real and positive energy variable e^. For 
the test system we assume, for simplicity, a single degree 
of freedom denoted by x which takes values x > 0. The 
energy of the test system is chosen to have the simple 
form £0 — H{x, A) ~ Ax, where A is the external forcing 
parameter. (Similar results were obtained for different 
energy functions, such as eo = H{x, A) = Ax^.) 

The parameter A is updated according to a set deter- 
ministic protocol A(t) in the time interval (0,r). During 
this time, the CS evolves with random sequential stochas- 
tic dynamics, with each site < i < L in the CS inter- 
acting with an arbitrary rate r with the site to its right 
by stochastically redistributing their total energy: 

4 = q{ei+ei+i) and e-+i = {l-q){ei+ei+i) . (13) 

Here < q < 1 is drawn from a uniform distribution 
and the prime denotes the values after redistribution [2^ . 
The update rule for i = (the test system) and the role 
of parameter A in the dynamics are explained in more 
detail as follows. Since the energy of the test system eq = 
Ax, the maximum possible value of x in the next update 
step is Xmax = {^0 + ^i) / X. Now we generate a random 
number Xr which lies uniformly between to Xmax and 
the update rules for the sites and 1 are: eg = Ax^ and 
e'l = (eo + ei — Ax^) where eg and e[ are the energies of 
the respective sites after redistribution. In this way we 
ensure that the transition to any of the states, arising due 
to the energy exchange, arc equally likely. One can easily 
check that the updating of sites i = and i ~ 1 agrees 
with Eq. [13] for this particular model. In accordance 
with standard definitions Q the work done on the test 
system is given by 

Large r compared to dX/ dt, the rate of change of external 
parameter, means a quickly thcrmalizing reservoir; in the 
limit 7' ^ oo for any finite dX/dt, we expect the reservoir 
to lose all memory and become Markovian. 

When X{t) = Ag = constant, we update the test sys- 
tem and the reservoir by repeatedly using the dynamics 
in Eq. [131 Clearly the total energy of the CS is con- 
served, and the dynamics satisfy detailed balance with 
respect to a measure that is uniform on a constant en- 
ergy surface, as required in a micro-canonical ensemble. 
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It is straightforward to show that in the thermodynamic 
Umit L — > oo, the equihbrium energy distribution of any 
site i ~ 0, . . . L takes the form 

Pie,) = /?e-^^' (15) 

where (3~^ = ^f=o^i/(^ + '^^'^ distribution of the 
CS is a product measure: Pcsii^j}) = Yli (^^~'^'^' ■ This 
allows us to calculate the partition function of the test 
system for a fixed value of A: Z{X) = e^^^^dx = 
(/3A)~^ with the free-energy given by F{\) = — /3^^ InZ. 

As the test system is driven, there is local heating in 
the region of its contact with the reservoir [sOl , with the 
consequence that the energy distributions of the reservoir 
sites near to the test system (e.g. i = 1, 2, 3 etc) deviate 
significantly from the equilibrium distribution as given 
in Eq. 1151 We study this numerically by employing the 
following driving protocol \{t): Starting at A(0) ~ Aq 
we increase the forcing parameter to A(t) = Aq + A A 
at intervals 5t in n\ = t/6t discrete steps of size 5\ = 
AX/nx. Between these updates the dynamics specified 
in Eq. [13] are carried out using standard Monte-Carlo 
methods. By discretizing Eq. [TH the work performed 
due to changing the parameter from Afc to Xk + SX at 
fc-th update step is defined as AWk = H{xk, Xk + SX) — 
H{xk, Xk) where Xk is the test system degree of freedom. 
The total work performed is = J2kti AW^. 

For fixed AA, a reversible process occurs in the limit 
0. For this to occur, the limit — > oo 
must be taken first, followed by either taking r — > oo 
or r ^ oo. All other protocols, in particular those with 
finite nx, are irreversible. In the following we compare 
systems by putting r ~ 1 (which sets the time scale) and 
varying the time interval r. 

The effect of local heating on the single-site energy 
distribution is shown in Fig. 31 where we plot the (in- 
tegrated) distributions P>(e) = P{e')de' for several 
sites for two values of r and with the same protocol for 
A(t). We see that far from the test system [i 3> 1) the 
reservoir is Boltzmann-distributed (as in Eq. [T5|l . but 
there is significant deviation near to the point of contact 
with the test system. As expected, for large r this effect 
decreases. In a sense, local heating is a measure of the 
non-Markovian nature of the reservoir as experienced by 
the test system. This notwithstanding. Fig. |4] demon- 
strates that even in the presence of local heating, the 
Crooks' theorem is clearly satisfied. The validity of the 
Crooks' theorem automatically implies that the Jarzyn- 
ski equality also holds. 

One of the chief uses of the Crooks' theorem and 
Jarzynski equality is to experimentally measure the (rel- 
ative) free energy of a test system. For this purpose, 
local heating is of central importance, since it has a pro- 
nounced effect on the number of experimental runs which 
must be carried out to accurately determine the free- 
energy difference AF = Fg — F4 [3l|. Our discussion 
will focus on the Jarzynski equality, but similar conclu- 
sions hold for the Crooks' theorem. In an experiment. 
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FIG. 4: Top panel: The integrated final [t — r) energy dis- 
tributions, -P>(e) for several sites, i for a short (r = 50) and 
long (r = 5000) experimental time. Bottom panel: The solid 
line is exp{(3W - PAF), where AF = \n[{\o + AA)/A], 
obtained analytically. We take r = 1 in an arbitrary unit and 
measure the time with respect to it. Results are obtained 
for L = 100, T = 1, Ao = 1, (5A = 1 and nx = 10, and by 
averaging over 10'' realizations. 



the quantity e ^ is measured once per run. The exper- 
imental average over Nj^ runs is thus 

, Nr 

" 1=1 

where Wi is the work performed in the /-th experimental 
run. Assuming the validity of the Jarzynski equality, this 
gives an estimate of the free energy difference measured 
from the experiment: 

PAFm = - In 6-/3^ . 

As we show, the more pronounced the local heating, the 
larger Nji must be in order for AFm to accurately esti- 
mate the actual free energy difference. To see this, note 
that the scale of the estimated error of the average of 
X = e-Z^CM'-AF) jg gj,^ by the width ax = y/iX"^} - {X)'^ 
of the probability distribution P{X) [s^l- The central 
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FIG. 5: Top panel: ax as a function of r, and Nr = 10*. Bot- 
tom panel: The fractional error in the measured free energy 
difference AFm, versus Nr/g\ where averaging is done over 
10^ realizations for each value of Nn/ax, and the straight 
line has slope —0.5. Results are for L = 100, T = 4, Ao = 2, 
AA = 1. 



limit theorem states that for Nji runs, the error in the 



estimate of (X) scales as axN^ for large Nr. It is 
straightforward to show that the error in the estimate of 
AF scales in the same fashion. Therefore, the smaller 
ax , the fewer the runs needed to obtain an accurate esti- 
mate of /S.F . This is especially relevant when the experi- 
mental protocol drives the system in a highly irreversible 
fashion. 

Fig. [5] shows the dependence of ax on t for two val- 
ues of n\. As is clearly seen in the data, the longer the 
time of the experiment, the smaller ax- In the absence 
of local heating, such as occurs in the t — > cx) limit, ax 
achieves a parameter-dependent constant value, as seen 
in Fig. [51 We note that since n\ is finite, the process 
is irreversible even for infinitely long experiments, with 
the consequence that ax does not reach zero. However, 
ax diminishes with increasing n\ , tending to zero in the 
reversible limit discussed above. Fig. [5] shows the devia- 
tion, AFm — AF, of the estimated free-energy difference 
from the actual value as a function of the expected scal- 
ing variable N[i/a\. Once again, we see that the larger 
the local heating, the more experimental runs are needed 
to evaluate the free-energy Note that by considering a 
system with maximal local-heating, an upper-bound on 
the number of realizations can be obtained. This can be 
realized by making sudden changes in the test system, 
thus effectively disconnecting the heat reservoir after an 
initial equilibrium distribution has been attained. 
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